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ABSTRACT 

The Plane Wave Approximation (PWA) is widely used in finite element 
analysis to implement the loading generated by an underwater shockwave. 

The method required to implement the PWA in NASTRAN is presented along with 
example problems. A theoretical background is provided and the 
limitations of the PWA are discussed. 


INTRODUCTION 


Background 

The finite element method is commonly used for analysis of structures 
exposed to underwater shockwaves. Modeling the structure using the finite 
element method is of less concern than loading the structure with a shockwave 
in a fluid medium. The fluid and the structure interact and the loading will 
be modified by the structure moving in the water. If the structure moves 
faster than the fluid can respond, then the density of the water near the 
structure will diminish considerably, causing, in simplified terms, a void in 
the water. This condition is known as cavitation and will cause the loading 
from the shockwave to be abated. As the structure slows, the fluid density 
will return to normal, or the void in the water will collapse due to ambient 
water pressure. Water closure, sometimes called a water hammer effect, occurs 
when the cavitation void closes and water slams together causing a momentum 
transfer. If this happens before the shockwave has passed the structure, 
then the shockwave loading will continue. 

In order to implement the underwater explosion shockwave loading, many 
finite element codes employ a load function known as the Plane Wave 
Approximation (PWA). The chief advantage of the PWA is that it is much less 
complex than other methods which employ modeling the fluid itself. This 
simplicity translates into faster computer run times with less memory 
requirements . 

Consider a water bounded node of a finite element model. Lets 
y = displacement of the node 
t = time 

= p 0 e' t/e be the incident shockwave 
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p r = (j> ( t ) be the reflected pulse off the model 

p a = mass per unit area associated with each node 
p 0 = density of water 

0 = time constant of pressure pulse (assumed exponential) 

A = water surface area associated with each node 
The equation of motion of the node is: 

P a 9 = P 0 e 8 + 4>(t) 


Continuity at the surface of the node stipulates that: 

_ _t 

Po°y = p o e 8 ■ 

Eliminating (J> using equations (1) and (2) results in: 

__t 

P a 9 + Po c Y ■ 2P 0 e 6 

The resulting force applied at the node is given by: 

r t 1 


Fit) = Ap a y 


2 PqG 


0 _ 


Po c Y 


( 1 ) 


( 2 ) 


(3) 


(A) 


The above equations were originally derived by Taylor in 1941 
(reference 1). The solution is exact as long as the water does not cavitate. 
The water may cavitate at the node or a short distance away from the node . 
Cavitation at the node will happen when: 

-■§ ( 5 ) 

p 0 cy k 2P 0 e 


The nodal force derived from the PWA can be stated by combining 
equations (4) and (5). 


F(t) 


2 P 0 e 6 - p 0 cy 


p 0 cy * 2 P 0 e 9 
p 0 cy < 2 P 0 e 9 


( 6 ) 


Objective 

The objective of this paper is to demonstrate the use of the PWA for 
underwater shockwave loading in COSMIC NASTRAN for transient analysis 
problems . 
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IMPLEMENTATION OF PLANE WAVE APPROXIMATION IN NASTRAN 

The PWA can be implemented in COSMIC NASTRAN without the use of the DMAP 
option. The procedure is facilitated by the use of scalar points and extra 
points, as these points will add neither mass nor stiffness to the finite 
element model . Consider a scalar dashpot which is implemented by the CDAMP4 
card. The dashpot is constrained to ground by leaving one of the coordinates 
on the CDAMP4 card blank. This will allow the dashpot to function 
independently of the structure which is being analyzed. The dashpot is given 
the value of pgc and loaded with an exponential shockwave: 

_ t 

Pit) = 2 P 0 e * (7) 

Equation (7) is best implemented with the TL0AD2 card. NASTRAN will not 
handle an instantaneous acceleration; therefore, equation (7) must be ramped 
initially to its peak value. The ramp is established with the TLOAD1 card 
coupled with the TABLED1 card. A ramp of 10 time steps to peak value was 
chosen using engineering judgment. The additional 10 time steps increased 
execution time only marginally while allowing sufficient ramp time to use 
the TL0AD2 card. Therefore, no attempt was made to fine tune the number of 
time steps in the ramp further. The TL0AD2 card is made to start on the 10th 
time step by using the DELAY card. The TL0AD1 and TL0AD2 card are executed 
together with the DLOAD card as shown in figure 1. From first principles, the 
velocity of the scalar point is: 


u ( t) 


scalar p t. 


2 Pge e 

P 0 C 


( 8 ) 


The wet node, extra point, and scalar point are equated using the TF 
card, as shown in figure 2. Each water bounded node must be assigned its own 
unique scalar dashpot, extra point and TF card. The TF relation defines the 
extra point velocity: 

pt. ~ scalar pt. & i *•) wot node 


The N0NLIN3 card is used to actually load the finite element model . 
Each wet node must have its own unique N0NLIN3 card and corresponding TF 
card. As shown in figure 2, the N0NLIN3 relation is: 


Fit) 


wet node 


P 0 CA[uit) ^ pt ] ; il ( t) ex pt > 0 

0 '• pt. * 0 


( 10 ) 


( 6 )). 


By using equations (8), (9), and (10) the PWA will result (see equation 


Fit) weC: n0( j e < 


^2P 0 e 6 - p 0 cyit)\ 


P 0 cyit)\ ; u(t) ex pt > o 

0 u< c hx. pt. * 0 


(ll) 
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ORIGINAL PAGE IS 
OF POOR QUALITY 


Schematic of PWA 



ILLUSTRATIVE PROBLEMS 


Piston Problem 

Consider the one dimensional problem of a rigid piston in a cylinder 
with a constant pressure on one side and an exponential shockwave in a flui 
medium on the other (figure 3). The water can cavitate at the piston or 
at some distance away from the surface of the piston. The piston problem was 
solved in closed form by Gordon and Handleton in 1985 (reference 2) and is 
categorized by the term: 

1 = E* ( 12 ) 

Po C0 

The piston can be modeled in NASTRAN with a single plate element 
constrained to move only in one direction. Four scalar dashpots are used for 
the four nodes of the plate element, and the dashpots are loaded with an 
exponential shockwave using the parameters shown in figure 3. The procedure 
described above is utilized to employ the PWA in NASTRAN. The results are 
shown in figure 4 in terms of displacement vs. time for varying values of X. 
NASTRAN was executed until the point of maximum displacement occurred. For 
X = 5.0 no cavitation occurs, and the PWA solution offers excellent agreement 
with the closed form solution of Gordon. 

When X = 0.5 cavitation occurs at the plate, the PWA differs from the 
closed form solution in terms of maximum displacement by approximately 36Z. 
This difference happens after 15 msec. Prior to this the PWA in NASTRAN has 
good agreement with the closed form solution. The difference could be 
attributed to the fact that the closed form solution has the additional 
impulse due to water closure. The PWA does not model the closure event. 

At X = 0.05 the PWA in NASTRAN and the closed form solution of Gordon 
differ by 156Z in terms of peak displacement. Cavitation occurs away from the 
surface of the piston, in the water itself. The PWA in NASTRAN does 
not indicate any cavitation as it only tracks the water edge at the piston. 

The PWA and the closed form begin to separate at approximately 4 msec, much 
earlier than when X - 0.5. Again, much of the difference could be attributed 
to the lack of a water hammer effect in the PWA. 

It is apparent that the accuracy of utilizing the PWA for the piston 
problem is dependent upon the value of X. Practically speaking, p 0 c remains 
fairly constant and 0 can only be varied over a small range. The term which 
can allow the most deviation is the mass of the piston p a . Thus, the PWA 
works well for a large mass term which does not allow cavitation. As the mass 
decreases, the effects of cavitation become more pronounced and the accuracy 
of the PWA decreases dramatically. 

Somewhere between X = 0.5 and 5.0 is the crossover value of X, where no 
difference exists between the solution of Gordon and the PWA. Assuming a 
linear relationship exists between X and the peak displacement difference 
between Gordon and PWA, the crossover can be found. Accordingly, for 
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Fig. 3. The Piston Problem 










X - 0.64 and greater, no difference in displacement should exist between 
PWA and the closed form solution of Gordon. 


Circular Flat Plate With Mass 

In this example a 2-inch thick flat circular plate is welded to a rigid 
annular backing structure. The backing structure constrains the outside 
of the plate but allows the inner section of the plate to move. In the center 
of the plate a dummy mass is welded which weighs approximately the same 
as the plate. The plate has air on one side, water on the other. A short 
distance from the plate is an explosive charge which will release an 
exponential shockwave when detonated. 

The experiment was conducted by the Underwater Explosions Research 
Division of the David Taylor Research Center and is shown schematically in 
figure 5. Velocity meters were placed on the mass to record the 
experimental response . 

A finite element model was formulated of the plate and mass using 
COSMIC NASTRAN. Using symmetry, only a quarter of the plate was modeled using 
plate elements. The mass was formulated using triangular plate elements and 
the plate with quadrilateral elements. The plate was modeled as pin connected 
at the edge. The mesh used is shown in figure 6. The PWA as described above 
was used to load the model. 

In this case the lowest value of X is 0.84. Although this problem is 
very different than the piston problem, if the problems were analogous we 
would expect no cavitation and the PWA to be an excellent method in which 
to employ the underwater shockwave loading. The results from the NASTRAN 
finite element model are plotted versus the experimental velocity in 
figure 6. The analytical solution offers good agreement with experimental 
results in terms of initial average acceleration, average deceleration and 
peak velocity. These results support the use of the PWA for this application. 


CONCLUDING REMARKS 

The PWA is a useful engineering method for analyzing the shock response 
of naval vessels from underwater explosions. The applicability of the method 
is dependent upon the presence of cavitation in the water. 

The PWA can be employed in NASTRAN for transient analysis problems 
without using the DMAP option. Standard bulk data deck cards can be used. 

The method requires the calculation for each wet node of several 
parameters (A,P 0 ,6) as well as an angular correction term if the shockwave 
and structure are not perpendicular. Additionally, several bulk data deck 
cards will be needed for each wet node. An accounting system must be set up 
in order to tie together the correct bulk data deck cards to the proper wet 
node. If the finite element model is of sufficient size then accomplishing 
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Fie. 5. - Test Arrangement 
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Fig*. 6^ - Comparison of Results from Experiment and P1JA 


the above can be quite a labor intensive task. A preprocessor program, which 
could scan the NASTRAN input deck and output the required bulk data deck 
cards necessary to implement the PWA, is recommended. 
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